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^ ' Abstract 
l> 

\^ I Survey for cosmology. It will determine the local mass function of Hi 

0^ ■ clouds, detecting several hundred per decade of mass. Each of these will 

. come with a redshift and, for the more massive clouds, an estimate of the 

Qh| velocity width. This will provide an ideal database for peculiar motion 

O ' studies and for measurements of biasing of galaxies relative to the under- 

lying matter distribution. 

a 

Surveys — Galaxies: luminosity function, mass function — large-scale struc- 
ture of the Universe 

x_ 

. In this paper, I wish to discuss the implications for cosmology of the Parkes 

Multibeam Southern-Sky Hi Survey (hereafter referred to as 'the Multibeam Sur- 
vey'). Most importantly, it will provide the first unbiased survey for Hi clouds, 
either associated with galaxies or truly isolated: I will show below that many 
thousand new objects should be discovered whose mass-function will place great 
constraints on models of galaxy formation. The second important feature of the 
Multibeam Survey is that it will map the velocity structure of the local Universe 
in great detail. When combined with accurate TuUy-Fisher distances, this will 
enable strong constraints to be placed on cosmological parameters, principally 
the density parameter. 

1 The mass function of Hi clouds 

I will start by estimating the baryonic mass-function of collapsed halos, assuming 
that baryons and dark matter are distributed equally over the sky {i.e. there is no 
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biasing). I will take parameters appropriate for the standard Cold Dark Matter 
(CDM) cosmology (others would give similar results): density parameter, Qq = 1; 
Hubble parameter, h = ifo/100kms~^Mpc~^ = 0.5; baryon fraction, = 0.064; 
normalisation, ag = 0.6. 

An analytic estimate for the number density of halos as a function of mass 
was first provided by Press & Schechter (1974). To estimate the proportion of the 
Universe which is contained in structures of mass m at redshift z, the density-field 
is first smoothed with a top-hat filter of radius R, where m = (47r/3)pi?^ and p 
is the mean density of the Universe. F{m, z) is then defined to be the fractional 
volume where the smoothed density exceeds some critical density 5c- Assuming 
a gaussian distribution, then 



where a is the root-mean-square fluctuation within the top-hat filter and erfc is 
the complementary error function. The key step was to realize that fluctuations 
on different mass- scales are not independent. In fact, to a first approximation 
Press & Schechter assumed that high-mass halos were entirely made up of lower- 
mass ones with no underdense matter mixed in. Then F must be regarded as 
a cumulative mass fraction and it can be differentiated to obtain the fraction of 
the universe contained in structures of a given mass. 



To convert this to a number density of halos per logarithmic mass-interval we 
simply multiply by p: dn/dlnm = p/. The main drawback of this approach is 
that, because of the above assumption of crowding together of low-mass halos 
into larger ones, it seems to undercount the number of objects. However modern 
techniques {e.g. Bond et al. 1992) give the same analytic form, simply scaled by a 
factor of two in normalization. The modified formula gives good agreement with 
numerical simulations {e.g. Lacey & Cole 1994) for 5c = 1.69 (as is appropriate 
for a spherical, top-hat collapse of density peaks. 

Figure ^ shows the predicted mass function for the baryonic content of ha- 
los and contrasts it with the observed stellar mass function. I have assumed a 
Schechter luminosity function. 



where C = 1.7 x lO^/i^ LoMpc-^ = lO^^h-^LQ, a = -1.25, and the stellar 
mass-luminosity ratio is mstar/L = 15h. 

It is apparent that the model predicts approximately the correct number den- 
sity of normal {L ~ L^) galaxies. However, it gives far too many halos at both 
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m / (h 'Mq) 

Figure 1: The mass function of halos in the standard CDM cosmology (sohd 
hne) compared with the observed luminosity function of galaxies (dashed line) 
assuming m/L = 15h Mq/Lq. Also shown is the Hi content of normal galaxies 
(dotted line). 

higher and lower masses. The reason for the former discrepancy was first ex- 
plained by Rees & Ostriker (1977) and by Silk (1977). They compared the ratio 
of the cooling time of the gas in proto-galactic halos to the dynamical time of the 
halo. In small objects the cooling time is always shorter than the dynamical time, 
thus gas can cool to form stars and hence a visible galaxy. In larger systems, how- 
ever, the cooling time exceeds the dynamical time. Then mergers may shock-heat 
the gas before it has time to cool, thus preventing significant star-formation. For 
this reason clusters of galaxies have little on-going star-formation except perhaps 
in a cooling flow deep within the cluster core. The dividing mass between these 
two regimes is highly sensitive to the gas fraction and cooling function but covers 
the range corresponding to the exponential cut-off in the luminosity function. 

The predicted excess of low-mass halos is harder to explain. It would seem 
that not all the cooled gas in proto-galaxies has formed visible stars. Until I 
researched this paper, it seemed possible to me that the missing gas resided 
inside the halos of low-mass galaxies in the form of Hi. The Hi content of normal 
spiral galaxies is insufficient (see, for example, the mass function shown by the 
dotted line in Figure |l| which is taken from the model in Briggs 1990) but I 
thought it possible that there may be a significant population of gas-rich dwarfs. 
However, the observations of the number counts of Hi clouds, described below, 
seem to rule this out. This is consistent with the optical luminosity function 
which, although it may miss many low-surface-brightness, predominantly low- 
mass galaxies is unlikely to have a faint-end slope as steep as the required value 
of a ~ —1.8. A more realistic explanation for the missing gas is that much 
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Figure 2: The detected density of Hi clouds {i.e. galaxies) of different mass, 
together with the expected sensitivity from the Multibeam Survey. For details 
see the text. 

of the proto-galactic interstellar medium was heated by an early generation of 
supernovae (and/or strong galactic winds) and expelled from the halo. If only 
a small fraction of Hi remains, however, or if it has fallen back into the galaxy, 
then it should be visible with the Multibeam Survey. 

Figure ^ shows the observed number density of Hi clouds in various surveys. 
The solid data points are taken from Briggs (1990). They all come from pointed 
observations towards different clusters, but including different proportions of fore- 
ground and background objects: solid squares, Leo group (Schneider et al. 1989); 
circles, Virgo (Hoffman et al. 1989); triangle, Hercules (Salpeter & Dickey 1985). 
Unfortunately these surveys are all highly biased. The only blind survey for Hi 
clouds which I know of is that of Kerr & Henning (1987) which has recently 
been analysed by Henning (1985). That gives a lower space-density as shown 
by the open squares in the figure. A realistic estimate of the true space-density 
is probably given by the lower locus of the points in the figure, some one to 
two orders-of-magnitude below the predicted curve if all the missing matter were 
in the form of Hi. As the observations must be shifted by about two orders-of- 
magnitude to the right to give good agreement with the predictions, this suggests 
that about one percent of the original baryonic mass of the halo persists in the 
form of Hi. 

The dashed line in Figure ^ shows the space density of Hi clouds at which 
the Multibeam Survey will detect 10 objects of a given mass. The sensitivity is 
almost three orders-of-magnitude higher than any previous survey thanks mainly 
due to the large volume of space which is covered. In making this prediction I 
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have have assumed that for a 5 sigma detection 
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mH.> 2.6x10^ Me Ut^J Lnn..^.-J ^ (4) 



lOMpc/ V 100 kms 

where the data have been binned to match the velocity width, Av, of the cloud. 
Supposing one per cent of the baryons to be in the form of Hi, then the TuUy- 
Fisher relation is 

/ \ 0.22 

Av ^ 440 I ^^^J^- I kms-\ (5) 

which gives a minimum density for the detection of clouds of 



^^-'■'^'°-'[-0)[mF^J ."-ft^Mpc- (6) 

The flat part of the sensitivity curve arises for Hi clouds which are detectable to 
the survey limit of 13 500 kms~^. 

On the basis of Figure I predict that the Multibeam Survey will detect 
several hundred Hi clouds per decade of mass. I would expect most of these to 
be associated with dwarf galaxies, but there may be a few truly intergalactic 
gas clouds. The mass and spatial distribution of these objects will be a major 
constraint on models of galaxy formation. 



2 Numerical simulations of Hi clouds 

The formation of bound objects by the growth and collapse of small density 
fluctuations in the early Universe is a highly complicated process. Although 
the Press-Schechter formalism, described above, gives a good estimate of the 
number density of objects, it tells us little about their spatial distribution or 
their internal structure. For this we have to rely on N-body simulations. Recently 
there has been a great improvement in the power of such simulations resulting 
primarily from two causes: firstly the introduction of massively parallel computers 
consisting of a large number of processors (typically a few hundred) each with 
their own memory and linked together by high-speed data channels, and secondly 
the development of sophisticated numerical algorithms able to take advantage of 
the new machines. Consequently pure N-body simulations {i.e. gravity only, or 
dark matter only) of a few tens of million particles and N-body, hydrodynamical 
simulations {i.e. a mixture of gas and dark matter) of a few million particles are 
now practicable. 

Studies of absorption lines in quasar spectra show that the high-redshift Uni- 
verse is full of Lya clouds, i.e. clouds of neutral hydrogen. Katz et al. (1995) 
simulate the production of these clouds in the standard CDM cosmology using 
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Figure 3: A histogram of the frequency distribution of Lya absorbing clouds 
per hne-of-sight {f{N) =d^A^/dA^Hid2;) in the numerical simulation of Katz et 
al. (1995). Error crosses and diagonal boxes show the observational constraints 
from damped Lya and Lyman limit systems, respectively. 

parameters similar to those described above. They use 64^ each of dark matter 
and gas particles with a gas particle mass of 7.75 x 10'^h~^ Mq. The comoving 
volume of the box is 11.11 /i~^Mpc and the box is evolved to a redshift of 2. 
A uniform photo-ionizing background is assumed to be present since a redshift 
of 6. Their paper contains a beautiful picture of the surface density of neutral 
hydrogen at the final time. It shows a dense network of interconnected filamen- 
tary structures studded with bright knots representing large concentrations of 
cold gas. They calculate the optical depth to Hi absorption as a function of 
frequency for a variety of lines-of-sight through the box. The resultant spectra 
are then analysed in a similar manner to quasar spectra to produce a histogram 
of absorption- line equivalent-widths. They are able to to reproduce the general 
form of the observed distribution over a wide range of equivalent widths, from 
10"^^-10^^cm~^. The low-equivalent- width systems arise from the filaments them- 
selves and from velocity caustics of the gas which is falling onto them. However, 
the high-equivalent- width systems, A^hi > 10^''cm~^, shown in Figure ^ occur 
when the line-of-sight passes through a lump of collapsed g galaxy (or 

more accurately a proto-galaxy as there is no star formation in the code). It 
can be seen from the figure that the model predicts too few absorption-line sys- 
tems. This may point to a deficiency in the standard CDM model, but is of little 
concern to us here. 

High-resolution simulations of this kind are very time-consuming and it is 
impractical to carry them forward beyond a redshift of 2 to the present day. Nor 
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would it be sensible to do so because of all the uncertainties in the physics of 
the intergalactic medium. However, the Multibeam Survey will anyway only be 
sensitive to column densities in excess of about lO^^cm"^ and we have just seen 
that these are associated with galaxies whose distribution can be determined 
in simulations of much lower resolution. This is one of the goals of the Virgo 
Consortium, a collaboration of mainly UK astronomers to carry out cosmological 
N-body hydrodynamical simulations of the formation of structure. We have been 
awarded time on the Cray T3D supercomputer at Edinburgh that consists of 512 
Dec alpha chips (of which 256 are usable at one time) each with 64 Mb of memory. 
We use the Hydra code developed by Couchman, Thomas & Pearce (1995) and 
available from |http: / /coho. astro. uwo.ca/pub/hydra/hydra.html| . Currently we 
are carrying out dark matter simulations with 17 million particles and dark matter 
plus gas simulations with 4 million particles. Only the latter are of relevance for 
this paper. 

Figure ^ shows the distribution of cold gas (T < lO^K) in one of our simula- 
tions a.t z = 0. It is a slice 7 /i~^Mpc thick through a box 70 h^^Mpc in width. 
The cosmological parameters are again similar to those given above, but the par- 
ticle mass is now 1.4 x 10^h~^ Mq. The box-size is well-matched to the size of the 
Multibeam Survey, although the mass-resolution is poorer than one would like 
and allows us to model the distribution of only the more massive galaxies. How- 
ever, I would expect an improvement of a factor of eight in mass-resolution within 
a couple of years. If one looks carefully then one can see numerous lumps of cold 
gas spaced along and at the intersection of filaments: these we associate with 
galaxies. It should be noted, however, that we have deliberately kept the physics 
in these simulations to a minimum. In particular, we have not attempted to form 
stars with all the associated feedback of energy into the interstellar medium. For 
this reason the Hi mass in the simulations is not representative of that in real 
galaxies. 

One of the main purposes of including gas in the simulations is to test the 
bias in the relative distributions of galaxies and dark matter. The former are 
more highly-correlated in space and are also moving more slowly than the latter. 
Moreover, the degree of biasing is dependent upon the mass of the galaxies, being 
larger for more massive systems. As an illustration of this. Figure ^ shows the 
relative distributions of moderate and high-mass galaxies in a test simulation. It 
will be very interesting to see from the Multibeam Survey whether there is a large 
population of dwarf galaxies filling the voids in the bright galaxy distribution. 

The degree of biasing is a major bug-bear of cosmology because it confuses 
the link between observations and theory: we see the galaxies, but the models 
predict only the overall distribution of matter. For example, one of the best ways 
to estimate the density parameter, Qq, is to measure the peculiar velocities of 
galaxies relative to the uniform Hubble expansion. The expected motions are 
proportional to and also to the overdensity of matter, that is 1/6 times the 
overdensity of light where b is the bias parameter. The Multibeam Survey will be 
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Figure 4: The projected distribution of cold gas at 2; = in a slice through a 
simulation of a critical-density, CDM universe (see text for details) . 
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Figure 5: Projections of the galaxy distribution in a test simulation: (a) m < 
0.5M„ (b) m > 0.5M,, where = 1.5 x W^h'^ Mq. 

an ideal database for peculiar motion studies. As well as detecting all large spiral 
galaxies out to more than 100/i"^Mpc, it will also measure their redshift and 
their Hi velocity width. When combined with infra-red photometry, this latter 
quantity will enable accurate TuUy-Fisher distances and hence peculiar velocities. 

Figure ^ shows 'wedge diagrams' for the real and velocity space distribution 
of galaxies drawn from an N-body simulation with Q = 0.3 and extent similar 
to the Multibeam Survey. The high density of galaxies which is expected in the 
Multibeam Survey provides an advantage over other sparser surveys: if one can 
find a well-defined void similar to those seen in the figure, then one can obtain a 
constraint on Qq which is independent of the bias parameter. This is because the 
underdensity in a void can never be greater than unity (whereas the overdensity 
in a cluster can in principle be anything). 

3 Conclusions 

1. The Parkes Multibeam Hi Southern Sky Survey will determine the mass 
function of Hi clouds in the local Universe. This will provide a strong 
constraint on models of galaxy formation. 

2. The spatial distribution of Hi clouds of differing mass will test ideas about 
biasing in the galaxy distribution. 

3. A complete survey of high- mass galaxies, plus TuUy-Fisher estimates of 
distances and hence peculiar velocities, will provide strong constraints on 
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Figure 6: 'Wedge diagrams' for the galaxy distribution in a simulation of spatial 
extent similar to the Multibeam Survey: (a) redshift space, (b) real space. 

the density parameter. 
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